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Abstract 

We are exploring the evolution of massive binary star populations for 
a nuclear starburst occuring in the conditions similar to the Milky Way 
galaxy, in its central few pc, on a time scale of 10 Myr. A computer 
code is applied allowing for computing, by using Monte Carlo simula- 
tions, the evolution of a large ensemble of binary systems, with proper 
accounting for spin evolution of magnetized neutron stars (NSs). Our 
results include the number of X-ray transients (NS + main sequence 
star), super-accreting black holes (BHs), and binaries consisting of a 
BH + supergiant, all as functions of time. 

We find that by 5 Myr after such a starburst one expects ~ 1 X- 
ray sources with a BH (Cyg X-l type), ~ 1 SS 433-like source and ~ 
54 transient sources with a NS, all to be within the central kpc. An 
interesting result that can be considered as a specific starburst feature, 
is that the ratio of the number of systems like SS 433 to the number of 
X-ray transients is about 1:100, compared to 1:1000, characteristic to 
the average ratio in the galactic field. The ratio of the number of X-ray 
sources containing a BH to the number of X-ray transients with NSs 
turns out to be a sensitive function of the age of the starburst, and its 
computed value ~ 0.03 is consistent with observations. 



1 Introduction 

X-ray observations of the Galactic center have revealed a number of 
energetic X-ray sources to be located in the innermost regions of the 
Galaxy (Churasov et al 1 ) ). Here, we apply Monte Carlo simulations 
to the evolution of a large ensemble of binary systems, with proper 
accounting for spin evolution of magnetized compact stars such as 
white dwarfs (WDs) and neutron stars (NSs) (Kornilov h Lipunov 

2 \ Lipunov 3 ), Lipunov et al 4 \ Lipunov 5 ), Lipunov et al. 6 ) ). In this 
poster, we will focus on the most prominent, from the observational 
point of view, representatives of the late stages of massive binary evo- 
lution, such as X-ray transients (NS in a highly eccentric orbit around 
a main sequence star, like A0535+26), super-accreting BHs (observa- 
tionally seen as SS 433 if in pair with a Roche lobe filling secondary 
component), and BH + supergiant binaries (like Cyg X-l). 

2 The model 

1. We consider only massive binaries, i.e. those having the mass of the 
primary (more massive) component in the range of 10 — 120 M . Only 
they can evolve off the main sequence during the time as short as 10 
Myr to yield compact remnants (NSs and BHs). 

2. Initial primary components' masses are distributed according to 
the Salpeter's mass function with an power law exponent a = 2.35. 

3. The initial mass ratio distribution was taken in a power law form 
f(q) oc q 2 . 

4. Initial eccentricity is assumed to be zero. 

5. After the collaps an additional kick velocity, Vkick, should be im- 
parted to the newborn NS. In the present calculations, the kick velocity 
was taken to be 75 km/s. 

6. We consider stars with a constant (solar) chemical composition. 

7. The process of mass transfer between the binary components is 
treated, when appropriate, as a conservative one, i.e. the total angular 
momentum of the binary system is considered to be constant. 

8. To get statistically significant results, the evoluiton of 300,000 bi- 
nary systems have been computed. Then we normalized the figures so 
as to be in agreement with the Tamblyn h Rieke 9 ) calculations of the 
number of massive OB-stars that survive ~ 7 Myr after the starburst. 



3 The results 



The evolution of some selected types of X-ray binaries during the first 
10 million years after the onset of star formation burst is presented in 
Figure 1. 

First we show (Fig. la) the number of X-ray transient sources con- 
sisting of a NS in an eccentric orbit around a massive secondary that 
acquired enough angular momentum during the mass exchange stage 
to become a rapidly rotating Be star. 

Then we show the evolution of a BH containing X-ray binaries of 
Cyg X-l type (Fig. lb) and of those with a superaccreting BH of SS 
433 type (Fig. lc). We note that the well known X-ray source IE 

1740-2942 (Churasov et al 1 ) ) may belong to SS 433-like binaries. 

One of the most interesting results concerns the relative number 
of SS 433-like binaries together with Cyg X-l-like binaries and X-ray 
transient sources. Our calculations yield for the ratio of these classes 
of sources to be ~ 1/50, which significantly exceeds this ratio for the 
entire Galaxy. 

The number ratio of a BH containing binaries (of both SS 433 and 
Cyg X-l type) to the X-ray transients with Be- stars is plotted in Fig. 
Id. It can be seen that this ratio is very sensitive to the time elapsed 
after the starburst and therefore it can serve as a tool to estimate the 
age of X-ray binaries at the Galactic center. 

In addition, the calculated a simple model spatial distribution of 
these sources (Fig. 2). 

4 Discussion and conclusions 

We compare our results with the observational data by using the 
GRAN AT X-ray observations of the Galactic center (Churazov et al 1 ) 
, Pavlinsky et al 8 ) ). 11 X-ray sources have been reported to be ob- 
served in the central region of the Galaxy (750 pc x 750 pc), including 
2 sources being classified as BH-candidates by their spectra, and 9 
being X-ray transients. 

As we noted above, the BH-candidates/X-ray transients ratio is a 
good indicator of the time passed after the onset of the starburst. The 
computed ration 0.2 at the age of 4 Myr and ~ 0.03 at the age of 5 Myr 
(which should be considered as a lower limit to the true value because 



we are not able to observe all the X-ray transients simultaneously) 
roughly corresponds to the observed ratio of ~ 0.2. Furthermore, in 
our model at the age of 7 Myr the X-ray transients occupy a region 
~ 2000 pc x 2000 pc in size (but more than the half of them are 
distributed in the region ~ 700 pc x 700 pc), as they acquire high 
velocities (after the faster evolving primary component explodes as a 
supernova) and move with an initially high velocity dispersion of about 
100 km/s. Hence, by considering the central 750 pc x 750 pc region, 
one will get the ratio closer to the observed value. 

To conclude, we note that the modelling of standard binary stellar 
evolution after a starburst at the Galactic center (Ozernoy 7 )) about 4-7 
million years ago yields the number of X-ray sources to be consistent 
with available X-ray observations of the Galactic center. 
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Table 1: Evolution of the selected types of X-ray sources 



Age 



2.0 



2.5 



3.0 



3.5 



4.0 



4.5 



5.0 



5.5 



6.0 



6.5 



7.0 



7.5 



8.0 



8.5 



9.0 



9.5 



BH/NS 



1.53E+02 



1.55E-00 



4.54E-01 



1.16E-01 



6.74E-02 



3.69E-02 



2.64E-02 



1.54E-02 



8.91E-03 



1.08E-02 



7.18E-03 



5.59E-03 



6.46E-03 



5.63E-03 



5.84E-03 



4.72E-03 



SS433 



2.87E+01 



3.02E+00 



1.25E+00 



8.25E-01 



6.84E-01 



5.79E-01 



5.61E-01 



4.91E-01 



3.33E-01 



4.04E-01 



4.56E-01 



3.16E-01 



3.33E-01 



3.51E-01 



3.86E-01 



3.51E-01 



Cyg X-l 



2.25E+01 



1.17E+01 



7.86E+00 



2.56E+00 



1.91E+00 



1.21E+00 



8.77E-01 



4.91E-01 



2.81E-01 



4.04E-01 



1.23E-01 



1.58E-01 



2.28E-01 



1.58E-01 



1.58E-01 



8.77E-02 



Transients 



3.33E-01 



9.47E+00 



2.01E+01 



2.91E+01 



3.85E+01 



4.85E+01 



5.45E+01 



6.39E+01 



6.89E+01 



7.45E+01 



8.06E+01 



8.47E+01 



8.68E+01 



9.04E+01 



9.31E+01 



9.29E+01 




Figure 1: The evolution of some selected types of X-ray binaries during the starburst at the 
Galactic center. 
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Figure 2: 2-D projection of the spatial distribution of some selected types of X-ray binaries 5 
Myr after the onset of the starburst at the Galactic center. Triangles- X-ray transients, square- 
SS433-like sorce, circle- Cyg X-l-like source. 





Initial distributions: 


9(a) 


~ 11 a, Abt & Levy (1976) 


9(M) 


~ M ' , Salpeter's law 


9(q) 


~ q\ oc q ~ 2 


9(1-0 


~ Kornilov & Lipunov (1983) 




Figure 3: Parameters of the model. 



